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Abstract 

We predict the differential cross sections for production of the X(3872) at the Tevatron and 
the Large Hadron Collider from both prompt QCD mechanisms and from decays of b hadrons. 
The prompt cross section is calculated using the NRQCD factorization formula. Simplifying as- 
sumptions are used to reduce the nonperturbative parameters to a single NRQCD matrix element 
that is determined from an estimate of the prompt cross section at the Tevatron. For X(3872) 
with transverse momenta greater than about 4 GeV, the predicted cross section is insensitive to 
the simplifying assumptions. We also discuss critically a recent analysis that concluded that the 
prompt production rate at the Tevatron is too large by orders of magnitude for the X(3872) to be 
a weakly-bound charm-meson molecule. We point out that if charm-meson rescattering is properly 
taken into account, the upper bound is increased by orders of magnitude and is compatible with 
the observed production rate at the Tevatron. 

PACS numbers: 14.40.Lb, 12.39.St, 13.85.Ni 

Keywords: Charm mesons, hadronic molecules, hadron colliders. 



I. INTRODUCTION 



The X(3872) is a cc meson that was discovered by the Belle Collaboration in 2003 through 
the decay B + — > K + + X [1]. Its existence was quickly confirmed by the CDF Collaboration 
through inclusive production of X(3872) in high energy pp collisions [2]. While some of the 
X(3872)'s at the Tevatron are produced by 6-hadron decays, most are produced promptly by 
QCD mechanisms [3]. A compelling case can be made that the X(3872) is a loosely-bound 
charm-meson molecule whose particle content is 

X = -^= [D*°D° + D°D*°) = (D*°D°) . (1) 
V2 

It has univeral properties that are determined by its binding energy only [4, 5] . They imply 
that the mean separation of its constituents could be larger than for ordinary hadrons by 
an order of magnitude or more. Given the identification of the X(3872) as a charm-meson 
molecule, an obvious question is whether the X(3872) could still be sufficiently robust to be 
observed at a hadron collider like the Tevatron. This issue has been brought into focus by a 
recent paper by Bignamini et al. [6] that derives an upper bound on the prompt cross section 
for a charm-meson molecule in terms of charm-meson-pair cross sections. They apply their 
upper bound to the X(3872) and conclude that the observed prompt production rate at the 
Tevatron exceeds their bound by orders of magnitude. 

The key to understanding how the X(3872) evades the upper bound of Ref. [6] is rescat- 
tering. In this paper, we will show that if rescattering of the charm mesons is properly taken 
into account, the upper bound of Ref. [6] is increased by orders of magnitude, which brings 
it into consistency with measurements at the Tevatron. Thus, despite its weak binding, the 
X(3872) is sufficiently robust to be studied at a hadron collider. Studies by the CDF Col- 
laboration at the Tevatron have provided the strongest constraints thus far on the quantum 
numbers of the X(3872) [7] and the most precise measurement of its binding energy [8]. At 
the Large Hadron Collider (LHC) at CERN, the data samples of the X(3872) will be much 
larger than at the Tevatron and thus will allow some of its properties to be studied with 
much greater precision. 

In this paper, we also predict the differential cross sections for production of X(3872) 
at the Tevatron and the LHC. We calculate the prompt production rates from QCD mech- 
anisms using the NRQCD factorization formalism [9], which expresses the cross section as 
the sum of parton cross sections for creating cc pairs with vanishing relative momentum 
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multiplied by phenomenological constants. We use simplifying assumptions to reduce the 
phenomenological constants to a single NRQCD matrix element, which is determined using 
an estimate of the prompt cross section at the Tevatron. We also calculate the rates for 
production of the X(3872) from 5-hadron decays using a method that gives the observed 
production rates of the J/ip from 6-hadron decays at the Tevatron. 

We begin in Section II by summarizing the case for the X(3872) as a loosely-bound charm- 
meson molecule and describing some of its universal properties. We also derive a factorization 
formula that determines the dependence of the production rate on the binding energy of the 
X(3872). In Section III, we use results from the CDF Collaboration to estimate the cross 
sections for production of the X(3872) at the Tevatron from both prompt QCD mechanisms 
and 6-hadron decays. In Section IV, we consider the relation between the production of the 
X(3872) and the production of pairs of charm mesons. We point out that if the effects of 
charm- meson rescattering are taken into account, the upper bound on the cross section for 
the X(3872) that was derived in Ref. [6] is increased by orders of magnitude. We also present 
an order-of-magnitude estimate of the X(3872) cross section and show that it is compatible 
with measurements at the Tevatron. In Section V, we use the NRQCD factorization formula 
for inclusive production to calculate the differential cross sections for prompt X(3872) at the 
Tevatron. We also calculate the differential cross sections for X(3872) from 6-hadron decays 
at the Tevatron. In Section VI, we calculate the differential cross sections for X(3872) at 
the LHC from both prompt QCD mechanisms and 6-hadron decays. Finally, we summarize 
our results in Section VII. 

II. THE X(3872) 

In this section, we summarize the case for the X(3872) as a loosely-bound charm-meson 
molecule whose particle content is given in Eq. (1). We also derive a factorization formula 
for production rates of the V(3872) that reveals how the rate depends on its binding energy. 

A. Universal Properties 

The only experimental information that is necessary to make the identification of the 
X(3872) as a loosely-bound charm- meson molecule is the determination of its quantum 
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numbers and the measurements of its mass. The quantum numbers of the X(3872) are 1 
which follows from 



• the observation of its decay into J/V'7, which implies that it is even under charge 
conjugation [10, 11], 

• analyses of the momentum distributions from its decay into J/ip-n + it~, which imply 
that its spin and parity are 1 + or 2~ [7, 12], 

• either the observation of its decays into D°D°7i° [13], which disfavors 2~ because of 
angular momentum suppression, or the observation of its decay into i/)(2S) 7 [14], 
which disfavors 2~ because of multipole suppression. 

The most recent measurements of the mass of the X(3872) in the <//?/> 7r + 7r~ decay mode 
[8, 15, 16] imply that its energy relative to the D*°D° threshold is 



The quantum numbers 1 ++ imply that the X(3872) has an S-wave coupling to D*°D°. Its 
tiny energy relative to the D*°D° threshold implies that it is a resonant coupling. Thus the 
X(3872) is an S-wave threshold resonance. 

A system of particles with short-range interactions that are fine-tuned so that pairs of par- 
ticles have an S-wave threshold resonance is predicted by nonrelativistic quantum mechanics 
to have universal behavior [17]. The universal properties of the system are determined by 
the pair scattering length a, which is large compared to the length scale 1/A set by the range 
of the interactions. The universal elastic scattering amplitude is 



where E is the energy relative to the scattering threshold and /x is the reduced mass of the 
pair. This scattering amplitude is applicable in the region \E\ < A 2 /(2/i), with errors that 
scale as (2/i|E|) 1 / 2 /A. If a > 0, one of the universal properties is the existence of a molecule 
with binding energy E x = l/(2/za 2 ). The universal wavefunction for this molecule is 



M x - (M B ,o + M D o) = -0.30 ± 0.40 MeV. 



(2) 



f(E) 



1 



(3) 



l/a + y/-2(j,E -ie' 



i>x(r) 



1 



exp(— r/a). 



(4) 
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The universal wavefunction in the momentum representation is 

~ Mr- (5) 

which is applicable for k < A. The wavefunction in the coordinate representation in Eq. (4) 
implies that the root-mean-square (rms) separation of the constituents is r x = a/y/2. 

To apply the universal properties of S-wave threshold resonances to the X(3872), we 
take the scattering length a to be that for charm mesons in the channel (D*°D°) + defined 
in Eq. (1). Assuming a > 0, the scattering length a = (2fiE x )^ 1 ^ 2 can be determined 
by using the binding energy from Eq. (2). Taking Ex = 0.30 ± 0.40 MeV as the input, 
we find that the charm mesons in the X(3872) have an astonishingly large rms separation: 
r x = 5.8l^° fm. A reasonable estimate for the momentum scale A associated with the range 
of the interactions between the charm mesons is the pion mass 135 MeV. The binding 

energy E x is very small compared to the corresponding energy scale m 2 r /{2jji) ~ 9.4 MeV. 



B. Factorization of the Production Rate 



Production rates of the X(3872) satisfy a factorization formula that is based on separating 
the momentum scale y / 2fiEx associated with the binding of the X from all the larger 
momentum scales of QCD [18], which include the pion mass and the charm quark mass 
m c . Our starting point for the derivation of the factorization formula is the differential cross 
section for a pair of charm mesons with energy near the threshold: 

da[D*°D°(k)} = fl^E/ d( f>n*D + .n \T[D*°D°(k) + all]| 2 ^^, (6) 

where k is the relative momentum of D*°D° in the center-of-momentum frame of the pair and 
"all" represents all the other possible particles in the final state. We have made a change 
of variables from the momenta of D*° and D° to k and the total momentum P of the 
charm-meson pair. The dependence on P has been suppressed in Eq. (6). In the remaining 
Lorentz-invariant phase-space measure rf0£>*D+aib D*D is treated as a single particle with 
mass M D *o + M D o m Mx- According to the Migdal- Watson theorem [19], the dramatic 
dependence of the T-matrix element in Eq. (6) on k due to final-state interactions resides 
in a multiplicative factor of the elastic scattering amplitude f(k 2 /2fi) given by Eq. (3). A 
factorization formula can be derived from Eq. (6) by dividing the T-matrix element by / and 
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then multiplying the right side by a compensating factor of \ f\ 2 . Since T/f is insensitive to 
the small relative momentum k, we can take the limit k — > in that factor. The resulting 
factorization formula is 

lM D->D« m = j |r [O -»B»(0) + al.]//(0)r x ,. 2 + W2 , )3V 

(7) 

The factorization formula for D°D*° is given by the identical expression. These factorization 
formulas hold for k < A, where A is the momentum scale set by the range of the interaction 
between the charm mesons. They indicate that the cross section integrated over k up to 
fc majf does not have the naive scaling behavior /c^ ax expected from phase space. Instead it 
scales like /c max in the range ^2/iE x < k nmx < A. 

In the corresponding factorization formula for the production of X(3872), the short- 
distance factor is the same as in Eq. (7). The long-distance factor can be deduced by using 
the fact that the spectral function for the resonance is proportional to the imaginary part 
of the scattering amplitude in Eq. (3): 

lmm = ^ (E + W)) + T 0tl_e { E). ( 8 ) 

The term with the factor 0(E) is the combined contribution from D*°D° and D°D*°. The 
term with the delta function is the contribution from X(3872). The long-distance factor 
in the factorization formula for D*°D° in Eq. (7) is the 0(E) term in Eq. (8) multiplied 
by dE/(4iT 2 ). In the corresponding factorization formula for X(3872), the long-distance 
factor is the 5-function term in Eq. (8) multiplied by dE/(2n 2 ) and integrated over E. The 
resulting factorization formula is 



a[X(3872)] = ^- £ f d<j> D * B+aU \T[D*°D°(0) + all]//(0)| 2 x (9) 

UX J Tlfl 

The dependence on the momentum P of the X(3872) has been suppressed, but P can be 
identified with the total momentum of D*°D° in the short-distance factor. 

1 /2 

Since the long-distance factor in Eq. (9) is proportional to E^ , the cross section goes 
to as the binding energy goes to 0. This is in accord with the common intuition that 
the production of a weakly-bound molecule should be suppressed in high-energy collisions. 
However the degree of suppression is often overestimated. The common intuition, as exem- 
plified by Ref. [6], is that the suppression factor should scale like E^ 2 . In the case of an 
S-wave threshold resonance, there is a much milder suppression factor proportional to E 1 ^ 2 . 
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A factorization formula equivalent to that in Eq. (9) was previously derived in Ref. [18]. 
In the factorization formula in Ref. [18], the long-distance factor is also proportional to E][ 2 . 
However the expression for the short-distance factor includes an explicit factor of A 2 , where 
A is the ultraviolet cutoff of an effective field theory for charm mesons with a large scattering 
length. That dependence on A is presumably cancelled by other terms in the short- distance 
factor that depend implicitly on A. Thus the expression for the short-distance factor in 
Ref. [18] can not be interpreted literally. The advantage of the factorization formula in 
Eq. (9) is that the short-distance factor is well-defined. 

III. PRODUCTION AT THE TEVATRON 

In this section, we use results from the CDF Collaboration to estimate the cross sections 
for X(3872) in pp collisions at the Tevatron from both prompt QCD mechanisms and from 
decays of b hadrons. 

Within a few months of the discovery of the X(3872) by the Belle Collaboration [1], 
its existence was confirmed by the CDF Collaboration through inclusive production in pp 
collisions at the Tevatron [2]. The X(3872) was also observed at the Tevatron by the .D0 
Collaboration [20], who showed that many of its production characteristics are similar to 
those of the ip{2S). The CDF Collaboration subsequently showed that for a specific data 
sample consisting of X(3872) in the decay channel J/tpn+n^ with a modest Pt cut, the 
production is dominated by prompt QCD mechanisms rather than 6-hadron decays [3] . The 
information given about this data sample can be used to estimate the cross section. 

In Ref. [3], the CDF collaboration observed both the X(3872) and the ip(2S) in a 
J/ip7T + 7T~ data sample with a cut Pt(J/iP) > 4 GeV as well as cuts on the pion momenta. 
For the X(3872), they applied an additional cut M(ir + 7T~) > 0.5 GeV to increase the 
signal-to-background ratio. For both the X(3872) and tp(2S), they measured the long- 
lived fraction / LL that come from the decay of b hadrons: f* L = 16.1 ± 4.9 ± 2.0% and 
/l L = 28.3 ± 1.0 ± 0.7%. The complementary fractions 1 — /ll are produced promptly 
by QCD mechanisms. Using the numbers of events reported in Ref. [3], we determine the 
product of the cross section for X(3872) and its branching fraction into J/ip7r + 7r~ for both 
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the prompt and 6-decay mechanisms: 

prompt [X (3872)] Br [X -> J/^ + ti~} = (0.0335 ± 0.0055) (e^/e x ) a tot ^(2S)}, (10a) 
^-dcca y [X(3872)] Br[X -> J/^tt+tt-] = (0.00643 ± 0.0023)(e^/e x ) <7 tota ih/>(2S)], (10b) 

where ex and are the efficiencies for observing X(3872) and i/)(2S) events in this decay 
sample. The ratio e^/e x is likely to deviate from 1 by tens of percents rather than factors 
of 2 [3]. The "total" cross section for ip(2S) in Eqs. (10) is the sum of the prompt and 
6-decay cross sections. The cross sections in Eqs. (10) should be interpreted as those for 
production of V(3872) and i>{2S) within the cuts used to define the two data samples. The 
cut pt{J/iP) > 4 GeV implies constraints on the transverse momenta of the X(3872) and 
ip(2S). The decay of X(3872) into J/-?/>7r + 7r~ proceeds through its decay into J/ipp*, where 
p* is a virtual p meson, with J/ip nearly at rest in the X rest frame. Thus the momentum 
components of X and J/ip differ only by the factor M x /Mj/^. We can therefore interpret 
the cross section on the left side of Eqs. (10a) as the prompt cross section for X — > J/%1)-k + h~ 
with pr(X) > 5.0 GeV. In the decay ip(2S) — > J/ipir + ir~, the momentum of the J/ip in the 
ip(2S) rest frame can range up to 0.5 GeV. Using Monte Carlo simulations, we have checked 
that the cross section for ip{2S) — > J /ip-n + n~ with p T (J/ip) > 4 GeV can be approximated 
by the cross section with p T (t/j(2S)) > 5 GeV within 12% accuracy. 

To obtain estimates of the cross sections for X(3872) from Eqs. (10), we need a measure- 
ment of the cross section for ip(2S). The sum of the prompt and 6-decay cross sections for 
ip(2S) in the phase-space region pr(tp(2S)) > 5 GeV and |y(^(2<S'))| < 0.6 has been mea- 
sured by the CDF Collaboration [21]. The product of this cross section and the branching 
fraction for ip(2S) into p + p~ is 0.69 ± 0.01 ± 0.06 nb. The branching fraction into muons is 
(7.5 ± 0.8) x 10~ 3 [22]. If we assume that the rapidity distribution for the X and the i/j(2S) 
are similar in the acceptance region of the CDF detector, we can insert this cross section 
measurement for vp(2S) into Eqs. (10). Using the estimate e^/ex ~ 1, we find 

<Vom P t[V(3872)] Br[X -> J/^tt+tt"] w 3.1 ± 0.7 nb, (11a) 
^-dcca y [V(3872)] Br[X -> J/^tt+tt-] w 0.59 ± 0.23 nb. (lib) 

We interpret these as the cross sections for X(3872) — > J/ipn + -K~ with pr(X) > 5 GeV and 
|y(X)| < 0.6. The error bars come only from the uncertainties in the i>(2S) cross sections 
and the numerical factors in Eqs. (10). They do not include the systematic error from setting 
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tip/tx = 1, which may be tens of percent. They also do not include the systematic error 
from approximating the cross section for tp(2S) — > J/ipir + iT~ with Pt(J/^>) > 4 GeV by 
the cross section with p T (ip(2S)) > 5 GeV. The estimate for the prompt cross section in 
Eq. (11a) agrees with the estimate given in Ref. [6]. 

Existing measurements provide some constraints on the branching fraction for the 
X(3872) into J/^tt+tt-. The Babar Collaboration has set a lower bound Br > 0.042 at 90% 
C.L. [23]. The sum of the measured branching ratios for J/-?/>7r + 7r~7r [10], [10, 24], 

D°D°n [15, 16] and ^(2S)-f [24] relative to J/ipn+ir- is 13.5 ± 2.9. If taken at face value, 
this implies an upper bound on the branching fraction into J/^it + it~: Br < 0.093 at 90% 
C.L. If we allow Br[X — > J/ip7i + ii^] to range from 0.042 to 0.093, our estimate of the prompt 
cross section for X(3872) from Eq. (11a) ranges from 72 nb to 33 nb, while our estimate of 
the 6-decay cross section from Eq. (lib) ranges from 14 nb to 6 nb. 

IV. CHARM-MESON-PAIR CROSS SECTIONS AND THE X(3872) 

In this section, we rederive the upper bound on the cross section for producing a charm- 
meson molecule presented in Ref. [6]. We point out that their calculation of the upper bound 
for the X(3872) is too small by orders of magnitude, because they do not take into acount 
charm-meson rescattering. We also present an alternative method for calculating the upper 
bound that is much more efficient than the method used in Ref. [6]. Finally, we derive an 
order-of-magnitude estimate of the cross section for producing X(3872) and show that this 
estimate is compatible with the observed prompt production rate at the Tevatron. 

A. Upper Bound on Prompt Cross Section 

In Ref. [6], Bignamini et al. derived an upper bound on the cross section for producing 
a charm-meson molecule in terms of the cross section for producing charm-meson pairs 
integrated over a region of small relative momentum. In the case of prompt production of 
X(3872) at the Tevatron, their upper bound was lower than the observed production rate by 
orders of magnitude, casting doubt on the identification of the X(3872) as a loosely-bound 
D*°D° molecule. 

The starting point in the derivation of the upper bound in Ref. [6] was an expression for 
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the inclusive cross section for the production of X(3872) in terms of its momentum-space 
wavefunction ip x (k): 



1 «^ f /-fcrnax j3t, 



2 



(12) 



The notation (D*°D°) + indicates a projection of the charm mesons onto the even-charge- 
conjugation component defined by the right side of Eq. (1). The integral over the relative 
momentum k has been restricted to the region \k\ < k nmx in which the integrand has 
significant support. By applying the Schwartz inequality to the integral over k in Eq. (12), 
one obtains an upper bound 

1 f , /" femax d 3 k 



1 r /-Kmax ML 

£ ita?/ ***** J WW, \n(D'°D° M k) + a»]f 

i'x(k) 2 . (13) 



all 

(2tt) 3 

The last factor is less than 1, because it is simply the incomplete normalization integral for 
the wavefunction. The remaining expression on the right side can be expressed as the sum 
of inclusive cross sections for producing pairs of charm mesons plus interference terms: 



isE/^y WW, \W° D °^ + all 



all 

- l<j[D*°D°(k < A; max )] + la[D°D*°(k < k max )] + (interference) . (14) 



In a high energy collision, an inclusive cross section is summed over many additional particles, 
which we have represented by "all" . The interference terms between the T-matrix elements 
for D*°D° and D°D*° should average to upon summing over all these additional particles. 
Thus the upper bound in Eq. (13) can be written in a very simple form: 

a[X] < ±a[D* D (k < A; max )] + \a[D°D*°(k < k max )]. (15) 

This inequality holds provided /c max is large enough to provide most of the region of support 
for the integral in Eq. (12). In the case of pp collisions, charge conjugation symmetry implies 
that the two cross sections on the right side of Eq. (15) are equal. 

In Ref. [6], the authors assumed that the region of support for the integral over k in 
Eq. (12) extends only up to the momentum scale set by the binding momentum (2/iEx) 1 ^ 2 
of the molecule. The central value of the binding energy from Eq. (2), E x = 0.30 MeV, 
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corresponds to a binding momentum of 24 MeV. The authors of Ref. [6] took the upper end- 
point of the region of support for the integral in Eq. (12) to be k max = 35 MeV. They used 
the Monte Carlo event generators Pythia and Herwig to calculate the distribution in the 
relative momentum k for prompt charm-meson pairs D*°D° at the Tevatron. The distribu- 
tion was normalized using CDF data on the production of D°D*~ at the Tevatron [28]. For 
the charm-meson-pair cross section integrated up to A; max = 35 MeV, they obtained 0.11 nb 
and 0.071 nb using Pythia and Herwig, respectively. These cross sections are more than 
two orders of magnitude smaller than our estimate of the prompt cross section for X(3872) 
from Eq. (11a), which is 33 — 72 nb. They concluded that the X(3872) is unlikely to be a 
charm- meson molecule. 

The flaw in the argument of Ref. [6] is that the authors ignored the effects of rescattering 
of the charm-meson pair. Rescattering can be important for relative momenta comparable 
to or smaller than the scale A set by the range of the interaction. For a generic molec- 
ular state whose binding momentum is of the order A, rescattering does not change the 
momentum scale. However an S-wave threshold resonance is characterized by a very large 
scattering length a that arises from a fine-tuned balance between the interaction strength of 
the constituents and the effects of rescattering. The rescattering effects are so strong that 
the cross section for elastic scattering saturates the unitarity bound for relative momentum 
in the range between (2fiE x ) 1 ^ 2 and A. These same effects allow charm mesons that are 
created with relative momenta as large as A to rescatter into small relative momenta of order 
\/2jlEx- For this reason, the region of support for the integral over the relative momentum 
in Eq. (12) extends up to the scale A. This dramatically increases the upper bound on the 
cross section given by Eq. (15). Now the charm-meson-pair cross sections calculated using 
Monte Carlo methods in Ref. [6] scale like fc^ ax from the phase space of the charm-meson 
pair. If we take m n as an estimate of the scale A and replace the value /c max = 35 MeV used 
in Ref. [6] by fc max = cA, the charm-meson-pair cross section of 0.07 — 0.11 nb is increased 
by a factor of 64c 3 . For reasonable choices of c, this upper bound can be larger than the 
estimated prompt cross section of 33 — 72 nb obtained in Section III. 

There is an alternative way to see that k max in the upper bound in Eq. (15) must be much 
larger than the binding momentum of the X(3872). By the Migdal- Watson theorem [19], 
the dramatic dependence of the T-matrix element in Eq. (12) on k resides in a multiplicative 
factor of the elastic scattering amplitude f(k 2 /2fx). If the region of support of the integral 
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over k in Eq. (12) was limited to the momentum scale (2fj,E x ) 1 ^ 2 , we could replace /(/c 2 /2/x) 
by the universal scattering amplitude in Eq. (3) and ipx(k) by the universal momentum-space 
wavefunction in Eq. (5). However the resulting integral over k in Eq. (12) is logarithmically 
ultraviolet divergent. This is a signal that the region of support for this integral is not 
limited to k of order ^2jiE X - It extends up to the momentum scale set by the range of 
interactions between the charm mesons. 

B. More efficient calculation of upper bound 

The upper bound on the cross section for the X(3872) in Eq. (15) is the cross section for 
charm-meson pairs integrated over the relative momentum k up to A; max . In Ref. [6], that 
upper bound was calculated for the case of prompt production of X(3872) at the Tevatron. 
We now present an alternative calculation of that charm-meson-pair cross section using a 
method that is much more efficient. In the analysis of Ref. [6] , the momentum distribution for 
charm-meson pairs was calculated by generating more than 5 x 10 10 2-to-2 parton-level events 
and then passing them through Pythia [25] or Herwig [26] for showering and hadronization. 
Only a tiny fraction of these events include a pair of nearly-collinear charm mesons with 
substantial pr- These specific events are generated by an outgoing gluon from the 2-to-2 
parton collision fragmenting into a charm-quark pair in the showering process, which in turn 
hadronizes most of the time into a pair of charm mesons. We calculate the charm-meson- 
pair cross section by generating gg — > gcc parton-level events with MadGaph [27] and then 
passing them through Pythia for showering and hadronization. This 2-to-3 parton-level 
process is the dominant mechanism for the production of a pair of charm mesons with small 
relative momentum and substantial Pt- In order to further increase the efficiency, only 
parton-level events for which the cc pair have relative momentum below 2 GeV have been 
generated. We checked that this reduction of the parton-level phase space does not affect the 
hadron-level cross section for charm mesons with small relative momentum k < 600 MeV. 
The main advantage of generating only gg — > gcc events is that the momentum distribution 
in the region of small relative momentum can be calculated with accuracy comparable to 
that in Ref. [6] by generating fewer events by about a factor of 10 4 . 

We follow Ref. [6] in normalizing the Monte Carlo distribution using measurements by 
the CDF Collaboration of the inclusive cross section for D°D*~ in pp collisions at the 
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FIG. 1: Cross section for the inclusive production of D°D*~ in pp collisions at y/s = 1.96 TeV. 
The cross section is integrated over bins of A<f> - the angle between the transverse momenta of D° 
and D*~ . The data points are measurements by the CDF Collaboration at the Tevatron [28] . The 
histogram is calculated using the Monte Carlo methods described in the text. 

Tevatron [28]. The CDF Collaboration measured this cross section in the phase-space region 
pr(D°), p T (D*~) > 5.5 GeV and \y(D°)\, \y(D*~)\ < 1. The cross section differential in 
- the angle between the transverse momenta of the two charm mesons - is shown in Fig. 1. 
In the analysis of Ref. [6], the Monte Carlo distribution for pairs of charm mesons was 
normalized to this data. Since the error bars are smallest in the bins of A0 closest to 90°, 
this means that in practice the Monte Carlo distribution in Ref. [6] was normalized to the 
data near 90°. In our analysis, as we are primarily interested in configurations in which 
two charm mesons are almost collinear, we normalize the Monte Carlo event rate instead 
to the measured rate in the first bin in A<p, which extends from to 15°. Compared to 
normalizing the rate near 90°, this increases the normalization by about 60% at the expense 
of a 20% error from the experimental uncertainty. Our Monte Carlo events for D°D*~ 
were generated via the parton-level process gg — > gcc, consistent with our procedure for 
generating D*°D° events. In order to increase the efficiency, we applied the parton-level 
cuts Pt(c), Pt(c) > 3.5 GeV, |j/(c)|, \y(c)\ < 2, and A0(c, c) < 90°, which do not affect the 
D°D*~ production rate in the region of interest. The rescaling factor - defined by the ratio 
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FIG. 2: Cross sections for the inclusive production of D*°D° in pp collisions at y/s = 1.96 TeV. (a) 
Differential cross section integrated over the region px(D*° D°) > 5 GeV and \y(D*°D°)\ < 0.6 as 
a function of the relative momentum k of D*°D°. The histogram is the hadron-level cross section 
calculated using the Monte Carlo method described in the text. The solid curve is the parton- 
level differential cross section integrated over the region pt(cc) > 5 GeV and |y(cc)| < 0.6 and 
normalized to the hadron-level cross section in the region k < 1 GeV. (b) Cross section obtained 
by integrating the parton- level differential cross section up to fc max . The dotted curves are simple 
phase-space distributions. 

of the experimental rate over the Monte Carlo rate in the first bin in A0 - is 1.62. This 
same rescaling factor was then applied to the distributions of D*°D° events. 

Having fixed the normalization of our Monte Carlo distributions, we proceed to calculate 
the production rate of D*°D° with small relative momentum at the Tevatron, run II. The 
differential cross section is displayed as a histogram in the relative momentum k in Fig. 2(a). 
The distribution, which was obtained by generating about 5 x 10 6 events, is as smooth as 
the distribution in Ref. [6], which was obtained by generating more than 5 x 10 10 events. 
For k less than about 300 MeV, the shape of the Monte Carlo distribution is consistent with 
k 2 , in accord with simple phase-space suppression. 

The results of our calculation using MadGraph and the event generator Pythia can be 
reproduced approximately using a simpler method that does not require any event generator. 
We assume that a D*°D° pair with small relative momentum k is formed from a cc pair with 
approximately the same relative momentum k. The resulting approximation for the D*°D° 
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cross section is proportional to the parton-level cross section for a cc pair: 

a[D*°D°(k < A; max )] w cr[cc(fc < fc max )] x F[cc -> D* ^ ]. (16) 

The factor F[cc — ► _D*°.D ] can be interpreted as the probability for a nearly-collinear charm- 
quark pair to evolve into a nearly-collinear D*°D° pair. The shape of this parton-level 
distribution agrees well with the hadron-level distribution shown as a histogram in Fig. 2(a). 
The normalization also agrees if we set F[cc — > D*°D°] = 0.142. The resulting distribution 
is shown as a solid line in Fig. 2(a). Alternatively, the normalizing factor can be determined 
from the CDF measurement of the cross section for D°D*~ in the first bin in A0 in Fig. 1. 
The resulting value is F[cc — > DD*\ = 0.118, which is smaller than the value obtained 
normalizing to the hadron-level distribution by about 17%. 

To obtain the charm-meson-pair cross section a[D°D*°(k < k mSLX )] that appears in the 
upper bound in Eq. (15), the k distribution in Fig. 2(a) must be integrated up to fc max . In 
Fig. 2(b), we show the cross section in Eq. (16) with F[cc — > D*°D°] = 0.142 as a function 
of /cmax- For the value /c max = 35 MeV used in Ref. [6], we obtain a cross section of 0.030 nb, 
which is smaller than that obtained in Ref. [6] using Pythia by a factor of 3.7. 

C. Estimate of the Prompt Cross Section 

The factorization formula for the production of X(3872) in Eq. (9) separates the short- 
distance effects associated with the production of the charm mesons from long-distance 
effects associated with their binding into the X(3872). The Monte Carlo methods used to 
calculate the upper bound on the prompt cross section at the Tevatron can be combined 
with this factorization formula to obtain an order-of-magnitude estimate of the prompt cross 
section. The Monte Carlo methods use an event generator, such as Pythia or Herwig, to 
model the hadronization process that produces the charm mesons. The tuning of the pa- 
rameters of the event generator implicitly takes into account typical hadronic cross sections. 
It does not take into account any unnaturally large cross section, such as that associated 
with an S-wave threshold resonance like the X(3872). 

The charm-meson pair D*°D° is expected to have a typical hadronic production rate as 
long as its relative momentum k is comparable to or larger than the momentum scale A set 
by the range of the interactions between charm mesons. According to the Migdal- Watson 
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theorem, the production rate for k < A is enhanced because of the factor f(k 2 /2fx) = (— 1/a— 
ik)' 1 in the T-matrix element. An event generator like Pythia contains no information about 
such an enhancement factor. When such an event generator is used, there is an implicit 
assumption that the T-matrix element T[(D*°D°) + (k) + all] does not vary dramatically 
for \k\ < A. If this assumption is applied to Eq. (6), we find that the "naive" cross section 
calculated using an event generator and integrated over \k\ < A can be approximated by 

a naivc [D*°D°(k < A)] « J- J2j # D .s +all \T[D*°D°(\k\ = A) + all] | 2 x (17) 

where the T-matrix element has been evaluated at \k\ = A. 

In the short-distance factor of the factorization formula for the X(3872) in Eq. (9), the 
term T/f is evaluated at k — 0. However this term is insensitive to the relative momentum 
for \k\ < A, so it can be approximated by its value when \k\ = A. The factorization formula 
can be expressed in the form 

(7[X(3872)] « I d <f>o*B +a n \T[D*°D°(\k\ = A) + all] A| 2 x (18) 

UX J 7T/i 

where we have used A 3> l/\a\ to approximate l//(A 2 /2/i) by A. The short-distance factor 
in Eq. (18) can now be eliminated in favor of the naive charm meson cross section in Eq. (17): 



<x[X(3872)] » a naive [D*°D°(k < A)] x 6?r v^^ , (19) 

If the naive cross section for producing charm mesons has been calculated using an event 
generator, this expression can be used to estimate the cross section for producing X(3872). 

In Section IV B, we used MadGraph and Pythia to calculate the naive cross section for 
D*°D° at the Tevatron as a function of /c max . Taking the central value Ex = 0.30 MeV for 
the binding energy from Eq. (2), our estimate of the prompt cross section from Eq. (19) is 
1.5, 5.9, and 23 nb for A = m n /2, m^, and 2m n , respectively. The upper end of this range 
of cross sections is compatible with the estimate 33 — 72 nb for the prompt cross section at 
the Tevatron obtained in Section III. 

We can use the estimate of the cross section for X(3872) in Eq. (19) to obtain a quantita- 
tive estimate of the upper endpoint k max of the region of support of the integral in Eq. (12). 
The upper bound in Eq. (15) can be approximated by the naive charm-meson-pair cross 
section <7 naivc [£>* /) (£; < k max )\ provided k max > A: 

a[X(3872)] < ( 7 naive [ J D*° J D°(A; < fc max )]. (20) 
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Demanding that the estimate in Eq. (19) is less than the upper bound in Eq. (20) even 
if y/2/iEx is as large as A, we obtain the condition k max > (67r) 1 / 3 A = 2.7 A. The naive 
cross section in Eq. (20) is shown as a function of /c max in Fig. 2. Setting /c max = 2.7 A, the 
upper bound on the cross section for X(3872) is 4.2, 32, and 230 nb for A = m n /2, m n , and 
2m T , respectively. This can be much larger than the prompt cross section observed at the 
Tevatron. 

V. NRQCD FACTORIZATION 

In this section, we discuss the NRQCD factorization formula for the inclusive production 
of X (3872), which expresses the cross section as a sum of parton cross sections multiplied by 
phenomenological constants. We use simplifying assumptions to reduce the phenonomeno- 
logical constants to a single constant that is determined from our estimate of the prompt 
cross section at the Tevatron. We then use the NRQCD factorization formula to predict 
the differential cross section for prompt production of X(3872) at the Tevatron. We also 
predict the differential cross section for production of X(3872) from 6-hadron decays at the 
Tevatron. 

A. NRQCD Factorization Formula 

Inclusive production rates for the X(3872) satisfy an NRQCD factorization formula, 
which separates the momentum scales of order m c and larger that are involved in the creation 
of the cc pair from all the smaller momentum scales of QCD [9]. The smaller momentum 
scales include the scale associated with the formation of the charm mesons and the scale 
\/2fiE x associated with their binding into the X(3872). The NRQCD factorization formula 
has the form [29] 

a[X(3872)]= £ a[cc n ] (<#), (21) 

n 

where the sum over n extends over the color and angular-momentum channels of a cc pair. 
The short-distance factors b are inclusive cross sections for producing a cc pair with negligible 
relative momentum in the channel n together with hard partons. They can be calculated 
using perturbative QCD. The dependence on the momentum P of the X(3872) has been 



17 



suppressed in Eq. (21), but P can be identified with the total momentum of the cc pair. The 
long-distance factor (O^) is proportional to the probability for the cc pair in the channel n 
to evolve into the cc meson H plus soft partons or hadrons. It is called an NRQCD matrix 
element, because it can be expressed as a matrix element of a local operator in an effective 
field theory for the cc sector of QCD called nonrelativistic QCD (NRQCD). The NRQCD 
matrix elements can be treated as phenomenological constants. They are universal, so once 
they have been determined by fitting experimental results, the NRQCD factorization formula 
can be used to predict the production rate in other experiments. 

The NRQCD factorization formula was originally developed for heavy quarkonium [29]. 
The relative sizes of the NRQCD matrix elements depend on the angular momentum quan- 
tum numbers of the quarkonium state, scaling as specific powers of the typical relative 
velocity of the heavy quark pair. The pattern of suppression is summarized by a rather 
intricate set of velocity-scaling rules. The NRQCD factorization formula can also be applied 
to pairs of charm mesons with small relative momentum [9]. In this case, there should be 
a hierarchy in the sizes of the NRQCD matrix elements according to their orbital angular 
momentum. The most important matrix elements should be the S-wave matrix elements. 
There are four independent S-wave matrix elements: (O* (* So)) , (0*( 3 Si)), (Og^So)), and 
(Og( 3 Si)). The subscript 1 or 8 on the operator indicates the color channel: color-singlet or 
color-octet. The argument is the spectroscopic notation 2S+1 Lj for the angular-momentum 
channel. Explicit expressions for these operators in terms of NRQCD fields are given in 
Ref. [29]. The S-wave matrix elements can be interpreted as probability densities for a cc 
pair created at a point in the specified state to form the hadron X multiplied by spin and 
color factors. For the four NRQCD matrix elements listed above, the products of the spin 
and color factors are 2N C = 6, 3(2iV c ) = 18, N% - 1 = 8 and 3(iV c 2 - 1) = 24, respectively. 
The large-scattering-length factorization formulas in Eqs. (7) and (9) imply that the rela- 
tive sizes of the NRQCD matrix elements for X(3872) should be the same as for D*°D°. 
Truncating the NRQCD factorization formula to the S-wave terms, we get 

<r[X(3872)] » afafSo)] (Of (%)> + a[cc 1 ( 3 S 1 )} (0?( 3 S 1 )) 

+a[cc 8 CSo)] (Of (%)> + %c 8 ( 3 Si)] (0*( ;i S 1 )). (22) 

A possible binding mechanism for the X(3872) is the accidental fine-tuning of the depth 
of the interaction potential between the charm mesons D*° and D° so that there is a bound 
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state very near the threshold. Since the mesons D*° and D° are color-singlets, the cc pair in 
the D*D system has equal probabilities 1/9 of being in a color-singlet state or in any of the 
8 color-octet states. It is therefore reasonable to assume that the probabilities for formation 
of the X(3872) from a color-singlet cc pair and from any of the color-octet cc pairs are equal. 
Given the normalizations of the NRQCD operators, this assumption translates into a ratio 
of 3/4 between color-singlet and color-octet NRQCD matrix elements: 

(Of( 3 ^)> « liOiCS,)), (23a) 
(O?CS )) « \iOfCSo)). (23b) 

We now consider the dependence of the NRQCD matrix elements on the spin quantum 
number, which can be spin-triplet or spin-singlet. In the hadronization stage of the formation 
of the X(3872), a cc pair evolves into the charm-meson pair (D*°D°) + with even charge 
conjugation defined in Eq. (1). Voloshin has pointed out that the cc pair in (D*°D°) + 
is necessarily in a spin-triplet state [30]. The approximate heavy-quark spin symmetry of 
QCD implies that transitions that change the spin state of heavy quarks are suppressed. 
The state (D*°D°) + is therefore much more likely to arise from a cc pair that is created at 
short distances in a spin-triplet state than in a spin-singlet state. As pointed out in Ref . [9] , 
this implies the suppression of the spin-singlet NRQCD matrix elements: 

{OfCSo)) « <0?( 3 Si)>, (24a) 
(O?CS )) « <<#( 3 Si)>. (24b) 



One possible mechanism for the binding of X(3872) that has not been excluded is an 
accidental fine-tuning of the mass of the charmonium state % c i(2P) to the D*°D° threshold. 
The cc state is transformed by its resonant coupling to D*°D° into a loosely-bound charm- 
meson molecule with the properties described in Section II. In this case, the relative sizes of 
the NRQCD matrix elements for the X(3872) may be governed by the velocity-scaling rules 
for 3 Pi charmonium states. The leading terms in the NRQCD factorization formula are the 
color-singlet 3 Pi and color-octet 3 Si terms [29]: 

a[X(3872)} » a[cc^Pi)\ (O? ( 3 Pi)> + a[cc 8 ( 3 ^)] (Of ( s Si)>. (25) 
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The color-singlet parton cross section cr[cci( 3 Pi)] and the color-octet NRQCD matrix element 
(0*( 3 Si)) both depend logarithmically on the NRQCD factorization scale A NR q CD in such 
a way that the dependence cancels between the two terms on the right side of Eq. (25). At 
any specific value of the momentum P of the cc pair, the color-singlet parton cross section 
cr[cci( 3 Pi)] can be made to vanish by adjusting Anrqcd- If the production is dominated 
by regions of P in which the ratio of o"[cci( 3 Pi)] to a[cc^Si)\ is roughly constant, then 
Anrqcd can be chosen so that the color-octet 3 S 1 term in the factorization formula in 
Eq. (25) dominates. Thus the S-wave factorization formula in Eq. (22) may be applicable 
even if the production of X(3872) is dominated by the formation of x cl (2P). 

B. Determination of the NRQCD matrix elements 

The NRQCD matrix elements in the NRQCD factorization formula in Eq. (22) are uni- 
versal. Once they have been determined by fitting experimental results, the NRQCD fac- 
torization formula can be used to predict the inclusive production rates of the X(3872) in 
other experiments. The estimate in Eq. (11a) for the prompt cross section at the Tevatron 
provides a single linear constraint on the products of the matrix elements and the branching 
fraction for X — > J '/ipn+Tr^ . In order to derive this constraint, we use MadOnia [31] to eval- 
uate the parton cross sections appearing in Eq. (22) at leading order in a s , which is order 
a 3 for pt(cc) > 0. We identify the 3-momentum P of the X(3872) with the 3-momentum of 
the cc pair in the center-of- momentum frame of the colliding hadrons. We set the mass of 
the charm quark to m c — 1.5 GeV, so a cc pair with vanishing relative momentum has mass 
3 GeV. The phase space region for the cc pair is Pt(cc) > 5 GeV and |y(cc)| < 0.6. We use 
the parton distribution function Cteq_611 [32]. The factorization and renormalization scales 
are set equal to the transverse mass (jp\ + 4m;?) 1 / 2 of the charm-quark pair. The resulting 
linear constraint on the NRQCD matrix elements from our estimate in Eq. (11a) for the 
prompt cross section at the Tevatron is 

Br[X -> J/^J7T + n-] {(0^( 3 S 1 )) + 0.159 (Of (%)) + 0.085 (Of (%)> 

+0.00024 (Of ( 3 Si))) = (2.7 ± 0.6) x 10~ 4 GeV 3 . (26) 

The error bar comes only from the statistical uncertainty in the estimate of the prompt cross 
section at the Tevatron in Eq. (11a). 
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Simplifying assumptions can be used to reduce the four independent S-wave NRQCD 
matrix elements in Eq. (22) to a smaller set. We will consider three different simplifying 
assumptions that reduce the four matrix elements to a single non-perturbative factor, which 
we choose to be (0^( 3 Si)). Our three simplifying assumptions are as follows: 

1. S-wave dominance . The X(3872) is equally likely to be formed from any cc pair that 
is created with small relative momentum in an S-wave state, regardless of the color or 
spin state of the cc pair. The assumptions on the NRQCD matrix elements are 

<0f( 3 Si)> = |<<#( 3 Si)>, (27a) 

(OfCSo)) = \(0^ 3 S l )), (27b) 

(O*eS )) = l(0^( 3 S 1 )). (27c) 

This is the same pattern of matrix elements as in the color-evaporation model 
for quarkonium production [33]. The linear constraint in Eq. (26) reduces to 
Br (Og ( 3 Si)) « 2.5 x KT 4 GeV 3 . 

2. Spin-triplet dominance. The X(3872) can be formed only from a cc pair that is created 
with small relative momentum in a spin-triplet S-wave state, but it is equally likely to 
be formed if the cc pair is in a color-singlet or color-octet state. The assumptions on 
the NRQCD matrix elements are 

<0?( s Si)> = liOfeS,)), (28a) 
<<WS>)> = (O*CS )) = 0. (28b) 
The linear constraint in Eq. (26) reduces to Br (Og ( 3 Si)) « 2.7 x KT 4 GeV 3 . 

3. Color-octet 3 S 1 dominance. The X(3872) can be formed only from a cc pair that is 
created with small relative momentum in a color-octet 3 Si state. The assumptions on 
the NRQCD matrix elements are 

(Of ( s 50> = (OfCSo)) = (OiCSo)) = 0. (29) 

This simplifying assumption was proposed by Braaten, who suggested that it might 
give a good approximation to the inclusive production rate of the X(3872) in 
high-energy hadron collisions [9]. The linear constraint in Eq. (26) reduces to 
Br (O* ( 3 S 1 )) w 2.7 x KT 4 GeV 3 . 
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FIG. 3: Cross sections for X(3872) — ► J/?/>7r + 7r~ in pp collisions at y/s = 1.96 TeV. The graphs 
are the transverse momentum (pr) distribution for \y\ < 0.6 (left panel) and the rapidity (y) 
distribution for px > 5 GeV (right panel). The curves are for prompt production assuming color- 
octet 3 <Si dominance (solid) or S-wave dominance (dotted) and for production from 6-hadron decay 
(dashed). 

Having obtained estimates of the NRQCD matrix elements for each of our three simpli- 
fying assumptions, we can now predict the differential cross sections for prompt production 
of X — > J/ipir + 7T~ at the Tevatron. The distributions of the transverse momentum and the 
rapidity are shown in Fig. 3. The solid and dotted curves correspond to color-octet 3 Si dom- 
inance and S-wave dominance, respectively. The curves for spin-triplet dominance cannot 
be distinguished from those for color-octet 3 Si dominance. For pt larger than about 4 GeV, 
our three simplifying assumptions on the NRQCD matrix elements lead to the same differ- 
ential cross section. The reason for this is that the order-a 3 parton cross section <3"[cc 8 ( 3 Si)] 
dominates over that for the other three S-wave channels at large Pt- In the color-octet 3 Si 
channel, the order-a 3 parton cross section includes a fragmentation contribution that de- 
creases asymptotically as dojdp\ ~ 1/px- For the other three channels, the order-a 3 parton 
cross section decreases more rapidly by a factor of or \jp%. They receive fragmentation 
contributions only at order a\ or oP s . Thus the differences between our three simplifying 
assumptions might be greater at large pt if higher order corrections in a s were included. 
These corrections are beyond the scope of this paper. As pt decreases below 4 GeV, the 
differential cross section for S-wave dominance continues to increase while that from the 
color-octet 3 Si dominance reaches a maximum and then begins to decrease. The reason the 
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S-wave dominance cross section increases is that the order-a^ parton cross sections in the 
color-singlet and color-octet 1 S channels diverge as pt — > 0. At very low transverse momen- 
tum, a fixed-order calculation for this channel is not reliable, and one has to resum higher 
order corrections in a s associated with soft-gluon radiation from the colliding partons. This 
resummation is beyond the scope of this paper. 

In the differential cross sections in Fig. 3, there is a 23% uncertainty in the normalization 
from the NRQCD matrix element. There are additional theoretical uncertainties associated 
with the NRQCD factorization formalism. The normalizations of the parton cross sections 
are sensitive to the value of m c and the renormalization and factorization scales. However 
the shapes of the differential cross section are much less sensitive. Thus variations in these 
parameters can be largely compensated by changes in the value of (0^( 3 Si)). 

C. Production from 6-hadron decay 

Beside prompt production, the other source of X(3872) in a hadron collider is the feed- 
down from decays of b hadrons. At the Tevatron, the rate from 6-hadron decays is smaller 
than the prompt rate but nevertheless significant. In the sample of X(3872) — > J/tpn + n^ 
events studied by the CDF Collaboration in Ref . [3] , the fraction of events from decays of b 
hadrons was 16.1 ± 4.9 ± 2.0%. 

To predict the differential cross section for X(3872) from 6-hadron decays at the Teva- 
tron, we follow a procedure similar to the one used in Ref. [34] to calculate the momentum 
distribution of J/ip from 6-hadron decays in high energy collisions. The differential cross 
section for X(3872) is the convolution of the cross section for producing a b quark, a frag- 
mention function for its hadronization into a b hadron, and the momentum distribution of 
the X(3872) from the decay of the b hadron. The Monte Carlo program MCFM [35] is used 
to generate the momentum spectrum of b quarks at next-to- leading order in a s . We set the 
mass of the b quark to m b = 4.75 GeV. We use the Cteq6_m parton distributions [32], with 
the factorization and renormalization scales set equal to the transverse mass {p\ + ml) 1 / 2 of 
the b quark. We use the Kartvelishvili fragmentation function [36] with exponent a — 29.1 
to describe the hadronization of the b quark into a b hadron. There are several ways in 
which we depart from the procedure used in Ref. [34]. We do not resum the logarithms of 
rrib/pT at large transverse momentum. To approximate the momentum distribution of the 
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X(3872) from the decay of the b hadron, we take the b hadron to be a B meson and we 
decay it into an X plus a kaon according to an isotropic distribution in the rest frame of the 
B meson. To normalize the differential cross section for X(3872) from 6-hadron decays, we 
use the estimated cross section for X(3872) — > J/ipiT + 7r~ at the Tevatron given in Eq. (lib). 
The normalizing factor can be interpreted as the product of the inclusive branching fraction 
for a b quark to decay into X(3872) and the branching fraction for X — > J /ip-n + it~: 

Br[6 -> X(3872) + any] Br[X -> J/^tt+tt"] = (1.9 ± 0.8) x 10~ 4 . (30) 

The error bar comes only from the statistical error in the estimate of the 6-decay cross 
section at the Tevatron in Eq. (lib). 

Having determined the normalizing factor in Eq. (30), we can predict the differential cross 
sections for producing X — > J/ipn + iT~ from b decays at the Tevatron. The predictions are 
shown as dashed curves in Fig. 3. The shapes of the distributions in pr and y are predicted 
to be similar to those for prompt production by the color-octet 3 Si mechanism except at 
low transverse momentum. The fraction of X(3872) from 6-hadron decay is predicted to 
increase from about 12% at pt = 5 GeV to about 37% at p T = 50 GeV. 

VI. PRODUCTION AT THE LHC 

In this section, we predict the differential cross sections for the production of the X(3872) 
at the Large Hadron Collider (LHC) from both prompt QCD mechanisms and the decays of 
b hadrons. We consider the production of X(3872) in the phase space regions of the ATLAS 
and CMS detectors and the LHCb detector. 

To predict the differential cross sections for production of X(3872) at the LHC, we use 
the same methods that were used to calculate the differential cross sections at the Tevatron 
shown in Fig. 3. We calculate the prompt cross section for X(3872) using the NRQCD factor- 
ization formula in Eq. (22). We consider the three simplifying assumptions for the NRQCD 
matrix elements introduced in Section VB: S-wave dominance defined by Eqs. (27), spin- 
triplet dominance defined by Eqs. (28), and color-octet 3 Si dominance defined by Eq. (29). 
The normalization of the prompt cross section is determined by the linear constraint in 
Eq. (26). We calculate the cc cross sections using the same parameters as in Section VB. 
We calculate the cross section for X(3872) from 6-hadron decays using the method described 
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FIG. 4: Cross sections for X(3872) — ► J/ipir + ir~ in pp collisions at y/s = 7 TeV. The graphs 
are the transverse momentum (pr) distribution for \y\ < 2.4 (left panel) and the rapidity (y) 
distribution for px > 5 GeV (right panel). The curves are for prompt production assuming color- 
octet 3 <Si dominance (solid) or S-wave dominance (dotted) and for production from 6-hadron decay 
(dashed). 

in Section VC. The normalization of the b decay cross section is determined by the product 
of branching fractions given in Eq. (30). 

We first consider the production of X(3872) in the phase-space region covered by the CMS 
and ATLAS detectors. In Fig. 4, we show the px distributions integrated over the rapidity 
range \y\ < 2.4 and the rapidity distributions integrated over the region p? > 5 GeV. 
The qualitative behavior of these curves is similar to those for the Tevatron in Fig. 3. 
The differences in the prompt cross section from different simplifying assumptions for the 
NRQCD matrix elements are negligible for px greater than about 5 GeV. The prompt and 
6-decay cross sections integrated over the region p T > 5 GeV and \y\ < 2.4 are predicted to 
be about 49 nb and 8.2 nb, respectively. The fraction of X(3872) events from 6-hadron decay 
is predicted to increase from 10% at pt = 5 GeV to 35% at pt = 50 GeV. However there 
are large uncertainties in the predictions. The statistical uncertainties in the normalizing 
factors are 23% for the prompt cross section and 39% for the 6-decay cross section. There are 
also additional theoretical errors from the masses of the charm and bottom quarks and from 
the choices of the factorization and renormalization scales. Although these uncertainties 
are large, they will partly cancel between the extraction of the normalizing factors from the 
Tevatron data and the calculation of the cross sections at the LHC. 
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FIG. 5: Cross sections for X — * J/ipir + ir~ in pp collisions at y/s = 7 TeV. The graphs are the 
transverse momentum (pr) distribution for 1.6 < ij < 5.3 (left panel) and the pseudorapidity (rj) 
distribution for pt > 0.5 GeV (right panel). The curves are for prompt production assuming 
color-octet 3 Si dominance (solid) or S-wave dominance (dotted) and for production from 6-hadron 
decay (dashed). 

We now move on to the production of X(3872) in the LHCb experiment, which is a 
forward detector in the pseudorapidity region 1.6 < rj < 5.3. In Fig 5, we show the pt 
distributions integrated over the pseudorapidity range 1.6 < rj < 5.3 and the pseudorapidity 
distributions integrated over the region > 0.5 GeV. The differences in the prompt cross 
section from different simplifying assumptions for the NRQCD matrix elements are small 
for pt greater than about 5 GeV, but they increase at smaller values of pr- The cross 
section for S-wave dominance is larger than that for color-octet 3 Si dominance by a factor 
of 2.5 at pt = 1.5 GeV and 6.5 at pt = 0.75 GeV. This leads to the large difference in 
the normalization of the rapidity distributions for these two cases in Fig. 5. The prompt 
cross section integrated over the region p T > 0.5 GeV and 1.6 < r] < 5.3 is predicted to be 
about 270 nb and 688 nb in the cases of color-octet 3 Si dominance and S-wave dominance, 
respectively. The 6-decay cross section integrated over the same region is predicted to be 
14 nb. Assuming color-octet 3 Si dominance, the fraction of X(3872) events from fe-hadron 
decay is predicted to increase from 2.2% at p T = 0.5 GeV to 9.0% at p T = 5 GeV and to 
20% at pt = 20 GeV. There are large uncertainties in these predictions, especially those 
that involve production at small pt- 
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VII. SUMMARY 



We have analyzed the production of the A (3872) at the Tevatron. We showed that the 
prompt production rate observed at the Tevatron is compatible with the identification of 
A(3872) as a loosely-bound charm meson molecule. In Ref. [6], the authors found that the 
production rate at the Tevatron was too large for a loosely-bound charm-meson molecule 
by orders of magnitude. The error in their analysis was that they did not take into account 
effects of rescattering of the charm mesons that are specific to an S-wave threshold resonance. 

We have predicted the differential cross sections for the A(3872) at the LHC from both 
prompt QCD mechanisms and the decays of b hadrons. The prompt cross section for A (3872) 
was calculated by using the NRQCD factorization formula. We used simplifying assumptions 
to reduce the nonperturbative factors to a single NRQCD matrix element (0*( 3 Si)), which 
was determined from an estimate of the prompt cross section at the Tevatron. The cross 
section for A (3872) from 6-hadron decay was calculated using a method that gives the 
correct differential cross section for J/ip from 6-hadron decay. The normalizing factor was 
determined from an estimate of the 6-decay cross section at the Tevatron. 

For the ATLAS and CMS detectors, we predict the production rate to be more than an 
order of magnitude larger than for the CDF detector at the Tevatron, given the same cut on 
the pt of the A and the same integrated luminosity. For the LHCb detector, where a much 
lower p T cut is possible, the production rate could be larger by another order of magnitude. 
Thus the LHC experiments should be able to collect very large data samples of A(3872). 

The large data samples of A (3872) at the LHC will allow precise measurements of various 
properties of this remarkable hadron. It should be possible to measure the branching ratios 
for decays of A(3872) with a good signature, including J/ipn + iT~, J/ipi, and ip{2S)^. It may 
also be possible to observe for the first time some rare decay modes of the A(3872), such as 
X c j(lP)7r + 7r~ [37, 38]. It should be possible to improve on the measurement of the binding 
energy of A (3872) by the CDF Collaboration [8]. It may even be possible to distinguish 
the universal line shape of the A(3872) [39, 40] from that of a conventional Breit-Wigner 
resonance. In conclusion, the experiments at the LHC are sure to add significantly to our 
understanding of the A (3872). 
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